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SUMMARY 

Wave guides play a significant role In microwave space communication sys- 
tems. The attenuation per unit length of the guide depends on Its construction 
and design frequency range. A finite element Galerkin formulation has been 
developed to study TM electromagnetic propagation In complex two-dimensional 
absorbing wave guides. The analysis models the electromagnetic absorptive 
characteristics of a general wave guide which could be used to determine wall 
losses or simulate resistive terminations fitted Into the ends of a guide. It 
Is believed that the general conclusions drawn by using this simpler two- 
dimensional geometry will be fundamentally the same for other geometries. 

In formulating the finite element solution, the continuous electromag- 
netic field Is divided Into a number of triangular discrete areas staked out 
by nodal points which are spread through out the field. The nodal points can 
be placed at any position desired. Thus, the finite element formulation can 
readily handle geometrical complications such as axial and transverse var1- 
at1ons*1n wall properties as well as cross-sectional area variations. The 
governing electromagnetic finite element equations and the appropriate bound- 
ary conditions are presented for a finite duct with known Incident modes. The 
entrance and exit boundary conditions are developed by coupling the finite 
element solutions to the eigen functions of an Infinite uniform perfect con- 
ducting duct. Example solutions are presented for electromagnetic propagation 
with perfect conducting and absorbing duct walls. 


INTRODUCTION 

Wave guides play a significant role In microwave space communication sys- 
tems. Depending on the material properties of the wave guide and Its design 
frequency range, the guide will attenuate a signal transmitted In the guide 
(ref. 1, p. 295). In general, losses In wave guides are undesirable; however. 
In power measurement applications, electromagnetic wave guides are sometimes 
fitted with fixed resistive terminations, movable vane or flop attenuators 
(ref. 1, p. 326) to absorb radiation and reduce reflections. To better Under- 
stand the microwave absorption process. In the present paper a finite element 
Galerkin formulation has been developed to study the power absorption charac- 
teristics of an electromagnetic wave guide. 

The numerical program to be presented herein Is complimented by Grant 
effort at Ohio State University, the University of Illinois, and Northwestern 
University. The work at Ohio State and Illinois Is primarily analytical In 
nature (ref. 2) while the research at Northwestern employs the transient finite 
difference technique (refs. 3 and 4). In addition. Rice (ref. 5) at NASA Lewis 
Is developing simple analytical correlation for electronic properties which 



produce the highest possible attenuation for a single mode In an Infinitely 
long duct. Such an analysis Is extremely useful for a numerical program, 
since the analysis gives a reasonable estimate of the electrical properties 
which will produce significant attenuation. 

The finite element solution developed herein models the variable property 
Helmholtz equation In a rectangular two-dimensional wave guide with complex 
structural wall boundary conditions and multiple mode Inputs. In the finite 
element solution, the continuous electromagnetic field Is divided Into a num- 
ber of triangular discrete areas staked out by nodal points which are spread 
throughout the field. Since the nodal points can be placed at any positioned 
desired, the finite element formulation can readily handle geometrical compli- 
cations such as axial variations In wall properties as well as cross-sectional 
area variations. This finite element solution bypasses the conventional 
eigenvalue problem with Its associated modes which have been considered In 
earlier works on electromagnetic propagation In ducts (ref. 6, p. 88). 

The Inlet and exit boundary conditions associated with the finite element 
analysis of electromagnetic propagation must be flexible enough to account for 
reflection of energy which can occur from wall changes Inside the guide or 
from abrupt exit terminations. The appropriate boundary conditions are simi- 
lar to those associated with flow gust disturbances, duct acoustics (ref. 7), 
and thermal wave propagation (ref./8). These problems occur In unsteady aero- 
dynamics. An addition modivatlon for the present study Is the continuing 
Interest at NASA Lewis In numerical modeling of unsteady aerodynamics and 
general wave propagation In ducts. 

In setting up the Inlet and exit boundary conditions, the variable prop- 
erty absorbing portion of the wave guide Is assumed to be joined by uniform. 
Infinitely long, and perfectly conducting entrance and exit ports. In the 
entrance port, a known source of electromagnetic energy Is assumed to be Inci- 
dent on the complex absorbing region which Is to be modeled by the finite ele- 
ment analysis. For a single dominant Input mode, Silvester, and Ferrar 
(ref. 6, p. 177) developed some mathematical constraints for the entrance and 
exit conditions. For this paper, however, the entrance and exit boundary con- 
ditions are developed for multimode propagation by coupling the finite element 
solutions to the complete set of eigen functions of a uniform conducting duct. 
In acoustics, numerical calculations using this boundary constraint were found 
to be In excellent agreement with experimental results (ref. 9). 

First, the governing electromagnetic equation and the appropriate bound- 
ary conditions are presented. Next, the finite element modeling Is presented. 
Immediately following the mathematical development, numerical validation exam- 
ples are presented and compared with known analytical results. Finally, the 
last set of example problems treat electromagnetic propagation with an absorb- 
ing boundary. 


NOMENCLATURE 

A/5 mode amplitude of plus going entrance waves, equation (29) 

Afj mode amplitude of reflected negative going entrance waves, 

equation (29) 
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magnetic flux density 

mode amplitude of plus going exit waves, equation (33) 

modified mode amplitude of plus going exit waves, equation (52) 

characteristic duct height 

dimensionless entrance height bg/b' 

dimensionless exit height b^/b' 

speed of light In vacuum 

electric field, equation (6) 

harmonic electric field vector, E(x,y,z) 

unit vectors In coordinate directions 

right hand side, equation (80) 

dimensionless frequency, equation (9) 

magnetic Intensity, equation (6) 

harmonic magnetic Intensity vector, H(x,y,z), equation (8) 

X component of magnetic Intensity 

I 

magnitude of Incident magnetic Intensity, equation (92). 

Incident magnetic Intensity, equation (91) 

reflected magnetic Intensity, equation (91) 

magnitude of reflected magnetic Intensity, equation (93) 

transmitted magnetic Intensity, equation (94) 

magnitude of transmitted mangetic Intensity, equation (94) 

maximum value of magnitude of magnetic Intensity In plotting 
domain, equation (106) 

minimum value of magnitude of magnetic Intensity In plotting domaV 
n, equation (106) 

normalizing magnitude of magnetic Intensity 

finite element approximation to Hx 

total number of line Increments, equation (46) 

Imaginary part 
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global stiffness matrix, equation (80) 

element stiffness matrix, equation (76) 

element stiffness coefficient, equation (77) 

wave number, equation (28) 

axial modal wave number, equation (30) 

element modal wave number, equation (78) 

dimensionless length, L'/b‘ 

position of material change, see equation (99) 

number of elements 

mode number, equation (48) 

number of nodes 

number of modes 

local interpolation shape functions 

mode number, equation (29) 

outward normal unit vector 

Poynting vectors, equation (87) 

total power transferred, equation (88) 

residual error, equation (60) 

real part 

length of line segment on boundary 
dimensionless time, equation (5) 
vector quantity, equation (64) 
weight, equation (58) 

dimensionless transverse distance, x'/b‘, equation (7) 
dimensionless transverse distance, y'/b', equation (7) 
dimensionless axial distance, z'/b‘, equation (7) 
damping coefficient, equation (96) 
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propagation coefficient, equation (96) 
reflection coefficient, equation (101) 
complex propagation constant, equation (92) 
dielectric constant, cVcq 
permittivity In entrance duct 
permittivity 
permittivity In vacuum 
conductance permittivity, equation (12) 
total permittivity, equation (12) 

Impedance, equation (13) 

permeability, w'/yQ 

dimensional permeability 

permeability In vacuum 

conductance, equation (7) 

transmission coefficients, equation (102) 

global unknown vector, equation (79) 

scalar, equation (62) 

angular velocity, equation (9) 

Subscripts 

scalar vector components 

Superscripts 

region 1 
region 2 

complex conjugation 
dimensional quantity 
transpose 
vector quantity 



GOVERNING EQUATION - AND MODEL 
Basic Model 

Consider an electromagnetic wave propagating from a transmitting antenna 
and Incident on a two-dimensional wave guide as shown In figure 1. As the wave 
encounters the Inlet, energy will be diffracted from the lip and redirected 
Into space, while some energy will couple to the duct modes and be transmitted 
Into the guide. Inside the guide, diffraction, absorption, and reflection 
will occur at variable property absorbing regions. Finally, energy will be 
transmitted down the duct and some will be reflected back out the Inlet. 

The present paper will focus on the Interaction of propagating duct modes 
with the wall as shown In figure 2. The model Is composed of three regions. 

In the uniform. Infinitely long and perfectly conducting entrance and exit 
regions, the exact solution of the governing differential equations can be 
easily written In terms of the duct modes. In the absorbing variable property 
region of figure 2, the finite element analysis Is employed to determine the 
field In the variable property nonuniform regions. As previously discussed, 
the transmitting antenna, (see fig. 1), sends electromagnetic waves Into the 
guide which are either reflected, absorbed, or transmitted In the absorbing 
region. Electromagnetic mode reflection at the Inlet to the duct absorbing 
region and the transmission at the outlet of the absorbing region are modeled 
by matching the finite element solution In the Interior of the nonuniformity 
to analytical eigenfunctions expansions In the uniform Inlet and outlet ducts. 
This permits a multimodal representation accounting for reflection and mode 
conversion by the nonuniformity (ref. 9). 


MAXWELL'S EQUATION 

The governing differential equations describing the propagation of elec- 
tromagnetic energy are the two Maxwell's equations. 

Vx E^ (x',y',z',t') = - ^ y'H^ (x',y',z',f) (1) 


V'xii^ (X' ,y' ,z' ,t' ) 


c*'E^(x',y',t') 


+ 


a 

at' 


4(x',y',z',t') 


( 2 ) 


In the foregoing equations, the prine, ', Is used to denote a dimensional 
quantity. These and all other symbols used In the report are defined In the 
Nomenclature. Vector equations (1) and (2) represent six scalar equations 
needed to determine the six unknowns scalars E^ (x' ,y ' ,z' ,t' ) , 


(x',y'.z',t 
and H 


). 


E^z(x- 


,(x',y',z',t'), H. (x',y',z',f). 


y'.z'.t'), H. 

^ 2 ^^' »y' subscript t inaicates time oependent fields, 

Later, when the assumption of harmonic time dependence Is assumed, the sub- 
script t will be dropped. 


t*- « <■ /, * '■ /* 

Indicates time dependent field; 


The nondImens lonal 1 zatlon begins with the dimensionless magnetic field 
H^(x' ,y' ,z' ,t' ) , permittivity c', permeability y', and electric field 

E|(x' ,y' ,z' ,t' ) and Introduces their nondimenslonal equivalents (no primes). 
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The scaling factor are the entrance duct height b', permittivity of a vacuum 

c', permeability of a vacuum v', velocity of light in a vacuum c', and 
0 0 0 
some magnitude of magnetic intensity Using these parameters, the 

dimensionless Maxwell's equations become 

vxE^(x,y,z,t) =• vR^(x,y,z,t) ( 3 ) 


7xH^(x.y,z,t) * a E^(x,y,z,t) 


+ ^ cE^(x,y,z,t) 


( 4 ) 


where c is the dielectric constant with the following dimensionless quanti- 
tative defined: 


t 




( 5 ) 


fl^(x,y,z,t.) 



E^(x,y,z,t) 


E'(x'.y.z'.t') 


( 6 ) 


a 


0 0 



( 7 ) 


As a further simplification, for harmonic steady state analysis it is 
usually assumed 

+jw't' - +1wt 

H^(x,y,z,t) » ^(x,y,z) e = H(x,y,z)e (8) 

where the dimensionless frequency is defined 

, O' “Lb* „ , 

f = , u = , w = 2* f (9) 

‘'0 ‘^o 


Similar assumptions for E are also employed. Substituting equations (8) and 
(9) into equations (3) and (4) yields 


9xE = -JwyH (10) 

7xH = jo)c.pE (11) 


where the total permittivity includes conduction 

CJ = C -Uc (12) 


The magnetic and electric fields are now considered to be harmonic spatially 
dependent variables. 
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NONHONOGENEOUS VARIABLE PROPERTY WAVE EQUATION 

The number of dependent variables can be reduced by combining Maxwell's 
equations (10) and (11) Into a single wave equation. This operation will lead 
to a 50 percent reduction of computer storage requirements. It Is also desir- 
able to develop a wave equation that could be used for varying media proper- 
ties so that no special treatment of the Interface between materials Is 
required. That Is, the same equations apply In the duct and In the absorber 
region and only the material properties are changed. 

Rewriting equation (11), 


VxH 

and taking the curl 



(14) 


(15) 


The constants ju are Independent of space and can be pulled out the curl 
operator In equation (15); however, cj niust remain Inside since cj now 
assumed a function of the spatial dimensions. 



Substituting equation (16) Into equation (10) yields our nonhomogenous 
governing wave equation. 


(16) 


vx 




(17) 


TM VARIABLE PR0PERTY20 WAVE EQUATION 

As customary In electromagnetics as well as fluid mechanics and heat 
transfer (ref. 10), simple Inlet profiles of the dependent variable are com- 
monly assumed. Herein, transverse magnetic (TH) waves will be assumed to 
represent the Input electromagnetic modes propagating down the entrance duct 
towards the nonuniform variable propertysectlon of the duct. The assumption 
Is made that only one component of the H vector will exist In the problem 
domain, that Is, 


H(y,z) 

where e^ Is a unit vector In the x 
fig. 3). The magnitude of the single 
only on the two spatial dimensions y 
would a z component of the magnetic 
property region. 


= Hj((y,z)ej( (18) 

direction (Into the paper as shown In 
X component of the vector H depends 
and z. The question remains, however, 
field be generated In the the variable 
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The validity of assumption (18) can be verified by consideration of the 
magnetic source equation (divergence equation); 

V • B » 0 (19) 

where B Is the magnetic flux density, a vector quantity defined by 

Bx = v(y.z)Hx (20) 

Recall, the above equation can be derived directly from Maxwell's first two 
equations (ref. 7, p. 280). Combining equations (19) and (20) and using the 
usual vector Identities (ref. 7, p. 578) equation 19 can be rewritten as 

v(y,z)v • H + H • 7v(y,z) = 0 (21) 

Substituting In equation (18) Into equation (21) yields 

u^[Hx(y,z)]^ Hx(y.z) ex • [^1^ ey ^ ej =0 (22) 

Since all the terms on the left hand side of equation (22) are Identical to 
zero, equation (18) Is a valid solution to the problem. Independent of property 
and geometry changes. 

Using equation (18) to evaluate the term In brackets In equation (17) 
yields 


vxH _ 1 a”x 

®T az 



(23) 


Next, the complete expression for the left hand side of equation (17) Is deter- 
mined by taking the curl of equation (23), as follows. 


7x 



(24) 


Finally, substituting equation (24) for the left hand side of equation (17) 
yields the scalar wave equation. 



(25) 


( 26 ) 


Equation (26) represents the governing wave equation to be solved by 
finite element theory. 
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UNIFORM DUCT ANALYTICAL SOLUTIONS 


The analytical solutions of equation (26) for wave propagation In the 
uniform perfect conducting Inlet and exit ports will be employed to give the 
proper termination boundary conditions for the finite element region. The 
analytical solution for TM waves traveling between parallel plates Is standard 
textbook material (ref. 11, p. 458) and quite simple In form. Consequently, 
these solutions will now be presented before the complete discussion of bound- 
ary conditions associated with the finite element solution In the variable 
property region. 


In the entrance and exit regions, the medium properties are assumed con- 
stant and real (0 = 0) such that the wave equation (26) reduces to 

v2Hx + k2Hx = 0 (27) 


where the free medium wave number k equals 

k = «.^/vc (28) 

In the analytical entrance region, the eigen value solution of equation (27) 
yields (ref. 11, p. 458) for the coordinate system shown In figure 3. 


H 


xa 




zn- 


(29) 


A y . i |( 7 

For the e'* time dependence used here, equation (8), the A e '' zn term 

” _ 2 

represents a wave propagating In the positive z direction while the An e zn 
represent a wave moving In the negative z direction (exactly opposite would 
occur had a e-J“t time dependence been assumed, see reference 11, p. 308. 


The axial wave number k^p equation (29) Is 




k 


Equation (31) represents the condition where the mode Is cut-off (nonpropa- 
gating). The negative sign Is chosen to produce mode decay In the positive 
direction for the A+ modes and mode decay In the negative direction for the 
A- modes. Since the characteristic duct height b* will be set equal to the 
entrance duct height, b^ will have a. numerical value of unity. 

The modal expression represented by equation (29) has been transacted to 
a total of N modes of the Infinite number possible. Thus, a total of 
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N unknown modal amplitudes iC, . . . have been Introduced. 
m \ c fifn m 

constraint equations w111 be required to determine each of these unknown 

reflection coefficients. The equations used to difine these coefficients will 

be Introduced In the following section on boundary conditions. 


The gradient of Hx In the axial direction will play a major role In 
the finite element solution to be presented. The gradient Is written as 

cos ^ ^ *^ b y^e"^*^zn^ 



Finally, In the exit port, a similar solution of the following form exits 



and 



(33) 





zn 


z 


(34) 


Here, only positive going waves are assumed and the eigenfunctions and wave 
number are expressed In the height of the exit duct bjj. 


BOUNDARY CONDITIONS 

A variety of boundary condition will be used In the finite element solu- 
tion of equation (26) for the model problem which Is displayed In schematic 
form In figure 2. Each of the required conditions will now be briefly 
discussed. 


Input Condition 

The analysis assumes a given number of propagating (see eq. (30)) A^ 
modes. These modes effectively set the level of the magnetic field In the 
finite element region and can be viewed as the equivalent Dirchlet boundary 
conditions required for ellptic boundary value problem as defined by 
equation (26). 
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Perfect Conducting Wall Conditions 

At a perfect conducting wall, the tangential component of the electric 
field vector Is zero (ref. 11, eq. 7-52a or ref. 12, eq. 1.69) 


Thus, for example. 


Et = 0 


Ex - Et » 0 


Along any of the horizontal perfectly conducting surfaces of figure 2. 

Since, the Hx equation (26) has been program for solution, the relationship 
between Ez given by equation (35) and Hx at the boundary must be developed. 
In particular, the finite element analysis well require Information on the 
gradient of Hx at a perfect conducting boundary. 

From equation (14), the Ez component of Maxwell equation can be written 
as (ref. 13, eqs. (1) to (4)). 


, ran 3H T 

- — - — 
Jwcj L3X ay J 


or since Hy a 0 


E - . -L- ^ 

^z “ “ jwc- ay 


Since E^ Is zero at perfect conducting horizontal boundary 


Now, since the finite element analysis will require Information about 
V Hx at Its boundary, let us expand at the upper horizontal boundary. 


VHx . n a VH^ 


/an an 

• ®y * liF * iT * 


ez J • ey 


Since only the dot product ey • ey gives a contribution. 


’"x • " ■ "aj <*' 

Since the electric filed tangency condition for this problem requires aH /ay 
Is zero (eq. (39)) ^ 


VHj^ • n a 0 


at the upper horizontal surface. It Is easily shown that this condition 
applies to all perfect conducting surfaces no matter their spatial orientation. 


12 


Continuity at Inlet and Exit 




In general, the tangential component of an H field Is continuous across 
an Interface (ref. 12, eq. (1.61)) except at perfect conductor where a surface 
current exists (ref. 11, eqs. 7-52(b)). Thus, the boundary condition becomes 

Hit = «2t (43) 


or the tangential components are Hx 

Hxa * Hx (z = 0; 0 < y < ba) (44) 


where Hxa the modal representation of magnetic field In the analytical 
Inlet region given by equation (29) and Wx represent the finite element 
approximation for Hx at the Interface. The hat over 'Pfx Implies an approx- 
imate finite element numerical solution to be discussed In detail In a follow- 
ing section. 


At the Inlet Interface, shown by the dashed line In figure 2, the H^a 
In the analytical region given by equation (29) must match the magnetic field 
defined by the finite element nodal points along the boundary Interface. Many 
possible matching methods can by employed for this boundary condition, such as 
point collocation, least squares or weighted residuals. A weighted residual 
approach was used herein with the weighting function equal to the 
eigenfunctions 



(y) - 


H(y)Jcos\^ b^ I 


dy 


z » 0 


(N^ equations, m s 1,2,3, N^) 


(45) 


Equation (45) represents separate equations; one for each coefficient 
defined In equation (29). 


For a particular weighting function cos ((m - l)iry/ba), a very simple 
numerical approximation for the line Integral can be written as 



(y^) - H 




(m - l)iry 


1 


4y^ 


» 0 (z 


0 ) 


(N„, equation, m = 1,2,3, ... Nm) (46) 

where the Ayi Is defined as shown In figure 4. Substituting In the value of 
Hj(a from equation (29) yields, moving known A+ coefficients to right 
hand side. 
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( 47 ) 



(Nm equations, m » 1,2, . . . Nm) 


where the mth equation is defined by the first cosine term on both sides of 
equation (47). In slightly different form 


m 


z 

n=1 




L i=l 


i 


m 


E‘n 

n®l 


i=l 

... V 1 

(z = 0) 

(Nm separate equations, m = 1,2, . . . Nm) (48) 


Z (‘" bj’* h) ' 

1 »1 


The terms on the right side represent the known forcing functions to the 
equation since all the values of are assumed known. The left hand side 
contains the Nm unknown reflection coefficients and the unknown nodal values 
of H1 at the finite element grid. Recall the modal equation used to specify 
the oncoming and reflected field was truncated to contain Nm unknown reflec- 
tion coefficients A", A", Ag, . . . Equation (48) represents one of the 

Nfp scalar equations required to determine all A^ reflective coefficients. 

At the exit, equation (45) becomes 



(y) - H(y)] cos 


z - 


» 0 


L a 0 


(Nm equations, m » 1,2,3, .. . Nm) 


or in terms of the modes 


(49) 


N 

m 


E 



(Nm separate equations, m = 1,2, . . . Nm) 
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(50) 



In contrast to equation (48), equation (50) has only the Bn values as 
unknown representing the positive going waves. No negative waves are assumed 
in the exit. In addition, the exit occurs at z equals L; however, the B+ 
coefficients have been redefined as follows 


n n on 


( 51 ) 


or 


6 


+ 

n 



+Jk 

e 


zn 


L 


In this manner, the exponential terms do not appear in equation (50). 


(52) 


Gradient at Inlet and Exit 

The last boundary condition required concerns the gradient of at 
the inlet and outlet of the nonuniforro finite element region. In the finite 
element analysis to be developed in the next section, a contour integral term 
will be developed which will contain the term • n. At a perfect conduct- 
ing wall, it was shown that this term is Identical to zero. However, along 
the entrance and exit interfaces, this term can be related to the modal equa- 
tion in the analytical region. At the inlet 

/dH 3H \ 

• " ■VaT % * »T 'zi • (-'z) ■ 


dH 


az 


(53) 


while at the exit 




(54) 


At the inlet and exit planes, continuity of the tangential component of 
the electric field is required across the interface (ref. 12, eq. 1.60). 
Employing equation (11) to express the tangential electric field in terms of 
the magnetic field (ref. 13, eq. 7-4) yields 



dH 


xa 

dz 


(55) 


For the special cases where the properties remain constant across the inter- 
face, equations (53) and (54) simplify to the following forms: 

dH 

VHx« n = - (inlet) (56) 
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and 


VHx* n = + (exit) (57) 

In general, the gradient of the magnetic field 1s continuous except when there 
Is a change In permittivity. 


FINITE ELEMENT THEORY 

The finite element formulation of the electromagnetic wave equation 1s 
now generated by using the Galerkin method (ref. 14) to obtain an Integral 
form of the variable property wave equation over the whole (global) domain D 
shown In figure 5. 


System Discretization 

The continuous domain D Is first divided Into a number of discrete 
areas staked out by the nodal points as shown In figure 6. Although the nodal 
points are shown evenly spaced In figure 6, the advantage of finite element 
theory Is that It allows the placement of nodes at any position desired. 
Another attractive feature of the finite element theory Is the ability to 
charge element properties In an easy manner. 

Next, the continuous magnetic field component Hx(y,z) will be approxi- 
mated (curve fitted) In terms of the nodal potential values Hxi located at 
yi , Zi , as shown In figure 6. The subscript 1 refers to any nodal point In 
this domain. The curve fit Is required since the variable properly wave equa- 
tion will be Integrated over the domain D shown In figure 5 and the curve 
fit will allow us to determine the value of Hx Inside the element. This 
contrasts with the simplest form of finite difference theory (taylor series) 
which usually only determines the values of the dependent variable at the 
lumped nodal points. 


GLOBAL WEIGHTED RESIDUAL APPROACH 

In the classical weighted residual manner, the magnetic field Intensity 
component Hx(y,z) Is curve fitted by expanding In terms of all the unknown 
nodal values Hxi(yi,z) and a series of basis (shape) functions, such that 

N 

Hx(y.z) * E W^(y,z)Hx^ = [W] |Hx| (58) 

where the basis or weight functions W^ (y,z) characterizes the spatial 
dependance of Hx(y,z) In terms of Hxi which represents the unknown value 
of the magnetic field Intensity component at the 1 nodal point In the global 
region. In the global approximation, the weight w^ has the property of 
being unity at node 1 and Identical to zero at all other nodes. As before In 
equation (44) the hat over the Hx(y.z) Indicates that It Is the approximate 
numerical solution to Hx(y,z). The nodes are numbered 1,2, .. . N and the 
global vector (Hx) represents the scalar values of the unknown magnetic 
Intensity component Hx at each node, such that In matrix form 
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since Hx(y.z) 1s only an approximation to the true solution Hx(y,z), 
substitution of equation (58) Into the governing wave equation (26) and Inte- 
grating over the domain D shown In figure 5 will not be equal to zero, as 
required by the exact solution but leave a residual error which Is defined 
as R 



— yH + JvH dydz = R 

Cf X xj 


(60) 


In accordance with the method of weighted residuals, the assumed basis func- 
tion and the distribution of errors R are forced to be orthogonal 
(R1 = 0) within the region by letting 



— yfi + 1 dydz = n 

Gy* ^ ^ / 


(1=1,2,... N equations) 


(61) 


Thus, there are n separate equations (written In compact form); one equation 
for each of the N nodal unknowns. In a direct analogy to the finite 

difference weighted residual control-volume formulation (ref. 15, p. 30), each 
of the above N equations represents a higher order difference approximation 
at the nodal point 1 where has a value of unity. 

By making use of the vector Identity of a scalar and a vector V 

y • (<;»v) = «|ry • v v • y<ii (62) 

and by letting t equal weight and V equal to 1/cj equation (61) 
can be expanded to 

• W,— H - — yH • yw, + dydz = 0 

1 Cy X Cy x 1 1 xJ •' 

(1 = 1,2, . . . N equation) (63) 
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Finally, using the divergence theorem, of Gauss (ref. 16, p. 79, eq. 4.7(b)), 



V) dydz = 



n) ds 


(64) 


where n is the unit outward normal vector to the boundary curve s, the 
first term in equation (63) can be converted into a surface integral such that 



(1 » 1 ,2, . . . N equations) 


(65) 


In effect, the second order differential equation has been reduced to a 
first order equation allowing the use of the weak formulation of the finite 
element theory. Thus in the Galerkin finite element approximation to follow, 
simple class Cq shape functions can be employed to approximate W^. Across 
an element shown in figure 6 for example, the class Co functions are only 
continuous in the dependent variable Hx and are discontinuous in slope. If 
the second order wave equation (61) is treated directly, then Hermitian calls 
Cl functions are required to approximate Wi (ref. 7). In these class C] 
functions, both the variable and its slope are continuous across a boundary. 


FINITE ELEMENT APPROXIMATION 


Both the specification of the global weighting function Wi and the glo- 
bal integration over the whole domain D in figure 5 required by equation (65) 
are not practical. However, the integration can readily be performed by sub- 
dividing the domain into smaller elements Ag, defining the global shape func- 
tion Wi in terms of the nodes of an individual element, integrating over an 
individual element and summing all the elements together. 

Equation (65) is valid over the entire domain D shown in figure 5 or any 
subdomain Ag, as represented by the area of a small triangular element embed- 
ded in the region as depicted in figure 5. To begin the finite element aspect 
of the weighted residual method, the domain D is assumed to be divided into 
M elements defined by N nodes, see figure 6. In this case, equation (65) 
can be written as 



7H. 




(1=1,2,... N equations) 


( 66 ) 


Where the properties are now given a subscript e to indicate that they may 
vary from element to element. Note, although the area integration in 
equation (66) applies to an individual element, the line integral is still 
defined over the global surface area, and as such, can be treated independ- 
ently from the local area integrations. 
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LOCAL SHAPE FACTORS 


Each element Is defined by the nodes around Its perimeter, while Hxi 
are defined at these nodes. However, the value of H^i inside the element 
are required to perform the integration. To represent the variation inside 
the elements of the field variable and its derivatives, local interpola- 
tion shape functions Ni(y,z) for linear triangles are written in both scalar 
and matrix form as 

yy.z) = H^®^(y.z) = ^ (y*z)»x3^ 


3 

= E 

M 


N 


(e) u(e) 


J 


xj 


= [N] 


(e) 


l«x( 


(e) 


(67) 


where {Hx> is the vector of nodal values of Hx for a general element 
e with subscripts 1, 2, and 3 representing the nodal positions as shown in 
figure 5. The subscripts 1, 2, and 3 take on the actual nodal number of 
figure 6 when applied to a specific element. 


The form of the local shape matrix [N] (®) depends on the type of element 
used. For the linear triangular element employed herein, the known value of 
Nj ^ is simple in form and can really be found in nearly every test on finite 

finite elements (ref. 6). Like its global counterpart W^, the local N^®^ 

has the property of being unity at node j and zero at the other nodes in the 
triangle. A graphical description of N is shown in figure 3 of reference 14. 

Replacing H by H' ' in equation (66), the new governing finite element 
equations became 




(1=1,2,... N equations) 


(68) 


6ALERKIN APPROXIMATION 

The weight W^(y,z) is now approximated by multiple values of n 1®^ (see 
ref. 14, fig. 3). ^Adapting the Galerkin approximation to the more general 
weighted residual approach assumes that 

W^(y,z) = Nj®^y,z) (69) 

in all elements containing the 1 node. For all elements which do not con- 
tain the node i, the weight is assumed zero not only at all other 
nodes (as required by the general definition of W^) but also at all values 
of y and z in the elements which do not contain the node i. 


(ei 

Recognizing that N) ' is zero for all elements not having the unknown 
Hxi associated with a particular element, the finite element equation (68) 
can now be written in compact form as 
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VH. 






u(e) 


7H 


fi^ ds 


0 


(i = 1 ,2, . . . N equations) 


(70) 


where are the known shape functions (ref. 6, eq. (3.05)). In general 
most of the M elements In the summation will have. zero contribution, since 
only If the 1 mode Is part of the element will have a nonzero value. 

Likewise, the line Integral can only contribute to a particular 1 equation. 
If the node Is on the boundary. For all Interior nodes, the line Integral 
term In equation (70) will not contribute, since ' will have a zero value 
on the boundary. 


SURFACE INTEGRAL AND BOUNDARY CONDITIONS 


The surface Integral In equation (70) can be rewritten along the boundary 
In figure 5 as 


f( 


1 X 


. n\ ds . / N,('> ^ . S) ds 

Entrance 

/ "S'* /"S'’ (’i 

ij. * “ ' nr r* * ^ 


nj ds 


(71) 


As shown In equation (42), the gradient dotted Into the outward normal Is 
Identical to zero along a perfect electrical conductor (PEC); thus, the last 
surface Integral on the right hand side of equation (71) Is Identical to zero. 
In addition, at the entrance, the axial derivatives In the line Integral at 
the entrance and exit planes represent the value of slope just Inside the 
finite element domain. They can be related to the slopes In the analytical 
region through boundary condition equations (56) and (57). Therefore, 



JL ,H(e) 

®Te 




ds 


^ ®Te 
Entrance 



dy 

z=0 


A.. ^ ®Te 


dz 


EXIT 


dy » 0 
z=L 


(72) 


Thus, the finite element equation (70) can be rewritten as 


t '"S'’ • 

e=l •'•'Ae ' 


7H<'>- . / N<'> ^ ^ 

Entrance ® 


dy 

z»0 




/XI 

N — sT 
' ‘Te “ 


Exit 

. . . N equations) 


dy = 0 
z=L 


( 73 ) 


(1 = 1.2, 
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where the gradients in the finite element region In equation (73) at the 
entrance and exit planes can be related to gradients In the entrance and exit 
ducts, equations (32) and (34) through equation (55). 


FINITE ELEMENT EQUATIONS 

Expanding Interms of the weighting function equation (67), and 

neglecting the surface Integrals for simplicity in writing (exact for central 
elements), equation (73) becomes 


e=l 

(1 s 1,2, . . . N equations) (74) 

(e) 

where the (Hx) local magnetic field Intensity component column vector can 
be pulled outside the area Integration since It does not depend on the area 
coordinates. 

Equation (74) could In principle be evaluated for each 1 to yield a set 
of N simultaneous equations for the N unknown value of However, 

this approach Is not readily Implemented on a digital computer. Rather, each 
element Is treated Independently and the contributions of a single element to 
all N equations defined by equation (74) are determine simultaneously from 
the following 



y/i 


i(e) 


Te 


[VN] 


(e) 




[N] 


(e)i 


dydz 3 0 



[vN^®^]- [N] dydz * 0 

(J = 1.2,3) 


(75) 


In contrast to the 1 subscript of equation (74) for all nodes In the global 
domain, here the subscript j sums only the three nodes of a particular ele- 
ment. In compact standard finite element form, equation (75) can be rewritten 
as 


[K] 


(e) 



» 0 


where 


K(e) 

U 




t N 


(e) 

1z 


■r 





dydz = 0 


with 


(76) 


(77) 


k 


2 

e 


“^»'e ‘Te 


(78) 
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For triangular elements, the evaluation of equation (77) Is quite simple and 
presented In many texts (ref. 18, p. 149, eq. (8.42) for the first two terms 
and page 45 for the last term). The surface Integral Is evaluated In a simi- 
lar manner. Shortly, a discussion Is given on how equation (76) Is used to 
build up equation (74). 


MODAL COEFFICIENTS 

In addition to the N finite element global equations (eq. 74) required 
to determine 'HxI values, additional equations are required for each unknown 
value Afj and BJJ which describe In the analytical entrance and exit 
region to the finite element domain. The unknown coefficient Ap and BJj 
couple to the finite element domain through the surface Integral terms In 
equation (73). For convenience In the finite element solution, the unknown 
nodal coefficients Af^ and Bjlj are Incorporated Into the global unknown 
vector {Hxl such that 

{*} ^ = Fa^ A 2 . . . Hx 2’ * * * ”x1 • • • ”xN ®1®2 • ‘ • Vl 

L m m J 


GLOBAL FINITE ELEMENT EQUATIONS 

Along with equations (48) and (50) which define the modes In the analyti- 
cal region, using the standard finite element procedures (ref. 18) the elemen- 
tal equation (76) can be assembled Into equation (74) to yield the foTlowIng 
global set of simultaneous equations which can be solved to obtain the unknown 
modal Intensities and unknown coefficient Ap and BJ{: 

[K] {«} = {F} (80) 

Here, the column vector F contains boundary condition Information which 
Includes Input modal condition A+ (right hand side of equation (48). The 
global stiffness matrix [K] Is the sum of the known local stiffness matrix 
[K](e) defined by equation (77). The solution of equation (80) by a banded 
Gauss solver yields i. 

Next, the solution of equation (80) for a number of element discretiza- 
tion patterns, absorber materials and nodal Inputs will be considered. 


RESULTS AND COMPARISONS 

For theory and code validation, first the finite element solutions will 
be applied to a number of problems where exact solutions exist. Next, the 
finite element solutions will be applied to a number of problem with various 
types of absorbers mounted on the walls. 

In all the examples to follow, the properties of the entrance region are 
assumed to be unity 
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In addition, the length of the finite element domain and height of the 
entrance duct are held at unity. 

For the linear element considered herein, the minimum required elements 
In the axial and transverse direction can be calculated from a formula devel- 
oped to give reasonable accuracy (ref. 7, eq. 48) using finite difference 
theory. 


Number axial element ^ 12 L g>/2« 


(82) 


Number transverse elements > 12 (n - 1) n > 2 


(83) 


where n Is the highest transverse mode to be resolved. In each example, a 
pictorial of the grid discretization Is shown directly above the plots of the 
spatial variations. 


Example 1. - Plane Wave Propagation In Perfectly Conducting Duct 

As the first example of the finite element technique, consider the 
problem of an Infinite long duct with perfectly conducting walls and with a 
plane wave propagating to the right from - Infinity with a frequency 
of 0 ) = 2ir. For this case, the wave continues propagating to + Infinity 
without producing any reflected waves. The exact analytical expression for 
Hx can be written (ref. 11, p. 457) as 

-Jk X 

|H^|= e = cos 2irx - 1 sin 2*x (84) 

As seen In figure 7, a comparison between the exact analytical results 
and the finite element analysis shows very good agreement between both analy- 
sis. The finite element analysis predicted A” coefficients of zero for the 
entrance reg1on^(no reflections) and unity for the b| coefficient while 
higher order B coefficients were appropriately near zero. 

Figure 8 applies to the same example only the magnitude of the Hx wave 
Is displayed where the magnitude Hx Is defined as 


|Hxl = yPe(Hx)2 + Im(Hx)2 


(85) 


Figure 9 displayed the axial flux of energy (Poynting vector) defined by 
(ref. 11 , eq. (8-69) 


Pave(y'Z) = J Re (E X H*) 


1 1 * '”x 

Pave 


( 86 ) 

(87) 


The total power radiated down the duct Is given by 


''t ■ 4 “ave 


( 88 ) 
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In all cases, the Inlet total power Is normalized as unity. Thus as expected 
In figure 9, the total power remains at a unity value throughout the duct. 

According to equation (87), the gradient of Hx Is required to determine 
the Poyning vector. In the analytical entrance and exit regions the gradient 
can, of course, be determined exactly from the known analytical expressions. 
However, In the finite element region some error Is Introduced since the gra- 
dient Is a constant throughout a triangular element (ref. 18, p. 30). This 
requires many small elements to approximate a rapid change In the variable 
Hx- However, If only the transmitted and reflected energy are required a 
coarse finite element grid will suffice, since the analytical expressions for 
the derivatives are known exactly. In figure 9, the small error In the finite 
element region (z * 0 to 1 ) results from the coarse elements, smaller elements 
will reduce this error. 


Example 2. - Third Propagating Mode In Perfectly Conducting Duct 


As the second example of the finite element technique, again consider the 
problem of an Infinitely long duct with perfectly conducting walls and the 
TH-j. (A*) wave propagating to the right from - Infinity with a frequency of 
u 5'^6.9 t 1 (f = 1.1). Again for this case the wave continues propagating to + 
Infinity without producing any reflected waves. The exact analytical expres- 
sion for TM 30 can be written according to equations (30) and (29) as 


k ^3 » 2ir X 1 


■1 V’ - (“)' = 


2.881 


-Jk z 

Hj^ a e = cos 2.881Z - J sin 2.881z 


(y - 0) 


(89) 

(90) 


where the frequency of the wave has been Increased over the value of 1 In 
example 1 to 1.1 so that the third mode Is now propagating while the fourth 
mode Is cut off. 


Again, as seen In figure 10 for the real and Imaginary parts of Hjj, 
and In figure 11 for the magnitude of Hx, the finite element and exact 
analytical results are In good agreement. 


Example 3. - Reflection, Transmission, and Absorption of Plane Waves 

As our next set of validation problems, consider the phenomenon of 
reflection which occurs when a uniform plane wave Is normally Incident on the 
boundary between regions composed of two different materials, as shown 
figure 12. The exact analytical results are given In standard textbooks. 
However, the exact analytical results have been slightly modified to account 
for the fact that the Incident plane wave contacts the Interface at z equal 
to L(j. 
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The exact analytical equations can be written In the following form 


"x” ■ "x 1 ‘ «xr 

(z < L^j) 

(91) 

"x 1 • 

(z < L^) 

(92) 

"xr- «ro 

(z < L^) 

(93) 

"x^’ ■ «xt ” "to ''^ 2 ' 

(z > L^) 

(94) 

where the magnetic field In region ( 1 ) 

Is composed of an 

Incident field 


H )(1 and a reflected field H^r whose magnitudes and propagation constants 
are defined In equations (92) and (93) with 

Y-| » + (c^^ = 0) (95) 

f 21 

In region (2), ' Is composed only of a transmitted magnetic Intensity 

deferred by equation (94) with 



(ref. 19, p. 388, eq. (32). For convenience In obtaining numerical results, 
Y 2 has been expressed In terms of 02 (real part) and 82 (Imaginary part). 
Or solving for a 2 and 62 (ref. 11, p. 364 P. 8-6 or ref. 12, p. 2-1) 



(97) 


(98) 


At z = L(|, the magnetic Intensity In region (1) must be equal to the 
magnetic Intensity In region (2); consequently, equations (91) and (92) are 
equated. Likewise, the Incident and reflected electric fields In region (1) 
must be equal to the transmitted electric field In region (2). The electric 
fields can be expressed In terms of the magnetic fields by employing Maxwell's 
equations. These two conservation equations can now be solved for the unknown 
reflected Hro and transmitted Hto wave magnitudes In terms of the known 
Incident magnitude 


£2 . r 

IO 


H 

H 


to 

I0 


(Y 
e 


^2^'"d 


(99) 

(100) 
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with the reflection r and transmission t coefficient defined as 


(ref. 13, eq. 5-97) 



2n, 


”2 * ^1 


(101) 


( 102 ) 


(ref. 13, eq. 5-98) 

The Impedance m and n2 have the usual convenient definitions (for 
positive values) 


n 


E 

-Jt 


(ref. 13, eq. 7-50) 


Jwc 


VxH 


T1 


^2 


E 

± 3£ 

1 VxE 



dz 


(103) 


(104) 


(105) 


The first case will consider a step change In permittivity from 1 to 4 
at z equal to 0.25. As shown In figure 13, the z equals 0.25 Interface 
occurs In the third row of elements which have been decreased In size to 
account for the smaller wave length In the higher permittivity material. As 
seen In figure 13, the finite element and exact analytical theories are In 
excellent agreement. Note the shift In slope of the real component at the 
Interface which Is automatically accounted for In the theory, because of vari- 
able property assumption In equation (26). Figure 14 shows similar good 
agreement between the finite element theory and the exact analytical results 
for the absolute magnitude of the magnetic Intensity The reflections 

from the Interface are clearly represented by the standing wave pattern ahead 
of the Interface. Observe that the magnitude of Increases In the mate- 
rial. The opposite effect will be seen In the next example. 

The second case considers a step change In permeability from 1 to 4 at 
z equals to 0.05. In this case, the Interface occurs Immediately after the 
very first row of elements shown In figure 15. As seen In figure 15, the 
finite element and exact analytical theories are In excellent agreement. For 
this particular example the slope Is continuous throughout both domains as 
expected since equation (55) Indicates that the slope Is continuous for 
changes In permeability. Figure 16 shows similar good agreement between the 
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finite element theory and the exact analytical results for the absolute magni- 
tude of the magnetic Intensity Again the reflections are clearly rep- 

resented by the standing wave pattern ahead of the Interface. In contrast to 
be previous example, the magnitude of decrease In the material. 

The final case considers a step change In conductivity from 0 to 2 at z 
equals 0.0. The jump In properties at z * 0 Is handled through equation (55) 
which Is employed In the entrance and exit Integrals of equation (73). As 
seen In figure 17 a slight difference In the magnitude between the finite ele- 
ment and analytical results occurs. The analytical solution was developed 
from a three region model (ref. 11, p. 347). A similar curve was obtained 
when the step change In conductivity was placed at z = 0.25 Inside the finite 
element region. Figure 18 shows a similar discrepancy for the Poynting vector. 
The error could not be eliminated with finer grid; therefore. It was assumed 
that the errors are a result of the discritizatlon process which must be 
resolved. 


Example 4. - Mall Absorbers 

Figure 19 shows a configuration with an absorber region along the upper 
wall. For a wall material with a permeability of 1 and total permittivity of 
l.-j2., the energy attenuation Is shown In figure 20. In figure 21, wall 
absorbers are placed at the top and bottom of the duct. For wall material 
with a permeability of 4.1 and permitivlty of l.-j 2.83, the wall absorption 
Is shown In figure 22. 


Figure 23 displays contour plots of the magnitude of the magnetic field 
for the absorber configuration shown In figure 21. The solid line drawing 
show the complete contours, while the symbol are plotted at selective points 
In the companion figure so that the actual magnitude of the contour can be 
readily determined. For convenience the magnitude of the magnetic field has 
been renormalized between 0 and 1 according to the formula 


l"x' - '%,n' 


contourl 


xmax 


|H 


xmln 


(106) 


where H^max the maximum value of magnitude of the magnetic field In the 
plotting domain and Hxmln Is the minimum value of Hx In the plotting 
domain. 


In figure 23(a) only one mode (plane wave) Is assumed In the modal expan- 
sion. In this case, 

|A*| = 1 .0 |A'| = 0.231 |B^| = 0.458 

and the value of the Poynting vector at the exit has fallen from 1 to 0.210. 

In figure 23(b) three modes are assumed In the modal expansion. In this case. 


= 1.0 

1^1 

= 0.1426 

Kl 

= 0.598 



-5 

1 +1 

1 -5 

= 0.0 

1^2 1 

= 0.475x10 

l«2 

1 = 0.98x10 
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jAgl = 0.0 |A"| » 0.664 IB 3 I = 1.123 

Because of the symmetrical nature of the magnetic field, the odd eigenfunc- 
tions have near zero coefficients. For this case, the Poynting vector at the 
exit has obtained a value of 0.358. 

In the contour plot example shown In figure 23(c), five modes are assumed 


In the modal expansion. In 

this case. 




Kl- 

1.0 

|A"| = 0.1246 1 


1 ■ 

0.6501 

lAjI - 

0.0 

lAgl = 0.186x10"^ 

|Bjl 


0.571x10 

|a;i ■ 

0.0 

|A~| = 0.6457 

|b;i 

= 

1.175 

ia;i • 

0.0 

|A~| = 0.376x10'^ 

GO 


0.398x10 

ia;i = 

0.0 

|A'| = 0.0145 


1 " 

0.0264 

For this case the 

Poynting 

vector at the exit was 

0 .- 

123. 



As seen In figure 23(c), the magnetic field has decreased to near zero at 
the walls but remains relatively high In the center of the duct. In effect, 
electromagnetic energy beams through the center of the duct. Future work will 
consider reducing the energy of this central beam. 


Finally, In the last contour plot example, seven modes are assumed In the 
modal expansion. This contour plot was for practical purposes Identical to 
the plot shown 1n figure 23(c). In this case, the modal coefficients are 


|a;i = ,.o I 

A‘| = 0.1204 

|a*| = 0.0 1 

A~| = 0.078x10 

IA 3 I = o-o 1 

A 3 I = 0.6416 

|aJ 1 = D-O 1 

|A"| = 0.555x10 

Kl = o-o 

|a'| = 0.0141 

|a;i = 0.0 1 

A" 1 = 0.799x10 

|a;i = 0.0 1 

A'l = 0.00286 


1b|| = 0.659 
IB 2 I = 1.185x10"^ 
[BjI = 1.160 
|bJ| = 0.951x10"^ 
IBjI = 0.0252 
|B^| = 0.959x10"^ 
|b^| = 0.00528 


and the normalized Poynting vector at the exit was 0.435, which Is very close 
to the value when five modes were assumed. In this seven mode expansion, the 
modal amplitudes of the first five modes are reasonably close to the values 
when a five mode expansion was assumed. In addition, the amplitude of the 
seventh mode Is very small. Therefore, a five mode expansion leads to 
convergence. 
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Example 5. - Thin Layer Absorbers 


The code has a zoom feature as shown 1n figure 24 for which the thin 
layered region Is expanded for ease In visualization. Now, any desired 
properly variation In the material can be conveniently reproduced. 

Figure 25 displays a multiple thin film region analyzed. The attenuation 
of the region Is shown In figure 26. 


CONCLUDING REMARKS 

The finite element method has been developed to handle the problem of 
electromagnetic propagation In a duct with varying wall properties and geome- 
tries. The numerical formulation Is relatively simple to use and appears to 
give very accurate results. For higher frequencies many wave lengths will 
exist In the duct requiring a large number of elements. With the expanded 
memory capability of the modern computer, the restriction does not represent a 
severe limitation of the method. 
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Figure 2.- Two dimensional wave guide finite element model. 













AXIAL COORDINATE, Z 


Figure 7. - A comparison of the axial magnetic intensity 

VARIATION ALONG THE LOWER WALL IN A UNIFORM DUCT WITH 
PERFECTLY CONDUCTING WALLS AS OBTAINED BY USING AN 
EXACT SOLUTION AND A FINITE ELEMENT SOLUTION FOR A 
PLANE WAVE (MODE 1) INCIDENT AT Z = 0 WITH U = 27T. 



AXIAL COORDINATE, z 

Figure 8.- A comparison of the magnitude of the axial 

MAGNETIC INTENSITY VARIATION ALONG THE LOWER WALL IN 
A UNIFORM DUCT WITH PERFECTLY CONDUCTING WALLS AS OB- 
TAINED BY USING AN EXACT SOLUTION AND A FINITE ELEMENT 
SOLUTION FOR A PLANE WAVE (MODE 1) INCIDENT AT Z = 0 
WITH U = 27T. 
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Figure 9. - A comparison of the axial flux of energy 
(POYNING vector) VARIATION IN A UNIFORM DUCT WITH 

perfectly conducting walls as obtained by using an 

EXACT SOLUTION AND A FINITE ELEMENT SOLUTION FOR A 
PLANE WAVE (MODE 1) INCIDENT AT Z = 0. 



AXIAL COORDINATE, Z 

Figure 10. - A comparison of the axial magnetic intensity 

VARIATION ALONG THE LOWER WALL IN A UNIFORM DUCT WITH 
PERFECTLY CONDUCTING WALLS AS OBTAINED BY USING AN 
EXACT SOLUTION AND A FINITE ELEMENT SOLUTION FOR MODE 3 
(TMjq) incident at Z = 0 WITH U = G.911. 
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AXIAL COORDINATE, Z 

Figure 11.- A comparison of the magnitude of the axial 

MAGNETIC INTENSITY VARIATION ALONG THE LOWER WALL IN 
A UNIFORM DUCT WITH PERFECTLY CONDUCTING WALLS AS OB- 
TAINED BY USING AN EXACT SOLUTION AND A FINITE ELEMENT 
SOLUTION FOR MODE-THREE INCIDENT (TMta) AT Z = 0 WITH 
U= 6.911. 



Figure 12.- A plane wave with normal incidence 

ON A PLANE BOUNDARY. 
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AXIAL COORDINATE, z 


Figure 13. - A comparison of the axial magnetic intensity 

VARIATION ALONG THE LOWER WALL IN A UNIFORM DUCT WITH 
PERFECTLY CONDUCTING WALLS AND A CHANGE IN PROPERTIES 
OF THE MEDIA AT Z OF 0.25 = 1,0 AND £2 = ‘•■0) AS 

OBTAINED BY USING AN EXACT SOLUTION AND A FINITE ELEMENT 
SOLUTION FOR A PLANE WAVE (MODE 1) INCIDENT AT Z = 0 
WITH U = 27T. 
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Figure 14. - A comparison of the magnitude of the axial 

MAGNETIC INTENSITY VARIATION ALONG THE LOWER WALL IN 
A UNIFORM DUCT WITH PERFECTLY CONDUCTING WALLS AND A 
CHANGE IN PROPERTIES OF THE MEDIA AT Z OF 0.25 (£^ = 

1.0 AND £2 = 4.0) AS OBTAINED BY USING AN EXACT SOLU- 
TION AND A FINITE SOLUTION FOR A PLANE WAVE (MODE-ONE) 
INCIDENT AT Z = 0 WITH U = 27T. 
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Figure 15.- A comparsion of the axial magnetic intensity 

VARIATION ALONG THE LOWER WALL IN A UNIFORM DUCT WITH 
PERFECTLY CONDUCTING WALLS AND A CHANGE IN PROPERTIES 
OF THE MEDIA AT Z OF 0.05 (Il.j = 1.0 AND P 2 ^ ‘•■0^ AS OB- 
TAINED BY USING AN EXACT SOLUTION AND A FINITE ELEMENT 
SOLUTION FOR A PLANE WAVE (MODE-ONE) INCIDENT AT Z = 0 
WITH U = 27T. 



Figure 16.- A comparison of the magnitude of the axial 

MAGNETIC INTENSITY VARIATIONS ALONG THE LOWER WALL IN 
A UNIFORM DUCT WITH PERFECTLY CONDUCTING WALLS AND A 
CHANGE IN PROPERTIES OF THE MEDIA AT Z OF 0.05 (|J^ = 

1.0 AND IJI 2 “ ‘*•0) AS OBTAINED BY USING AN EXACT SOLU- 
TION AND A FINITE ELEMENT SOLUTION FOR A PLANE WAVE 
(MODE-ONE) INCIDENT AT Z = 0 WITH U = 27T. 




AXIAL COORDINATE, Z 

Figure 17. - A comparison of the magnitude of the axial 

MAGNETIC INTENSITY VARIATIONS ALONG THE LOWER WALL IN 
A UNIFORM DUCT WITH PERFECTLY CONDUCTING WALLS AND A 
CHANGE IN PROPERTIES OF THE MEDIA AT Z OF 0.0 AS OB- 
TAINED BY USING AN EXACT SOLUTION AND A FINITE ELE- 
MENT SOLUTION FOR A PLANE WAVE (MODE-ONE) INCIDENT AT 
Z = 0 WITH U = 27T. 
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Figure 18 . - A comparison of the axial flux of energy 

(POYNTING VECTOR) VARIATION IN A UNIFORM DUCT WITH 
PERFECTLY CONDUCTING WALLS AND A CHANGE IN PROPER- 
TIES OF THE MEDIA AT Z OF 0.0 AS OBTAINED BY USING 
AN EXACT SOLUTION AND A FINITE ELEMENT SOLUTION FOR 
A PLANE WAVE (MODE-ONE) INCIDENT AT Z = 0 . 





AXIAL COORDINATE, Z 

Figure 19.- Discretization of air filled wave guide 

WITH ABSORBERS MOUNTED ALONG THE UPPER WALL. 
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AXIAL COORDINATE, Z 

Figure 20.- Effect of upper wall absorber on the magni- 
tude OF the axial flux of energy (Poynting vector) for 

A THREE-MODE MODAL EXPANSION IN THE ENTRANCE AND EXIT 
DUCTS. 
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AXIAL COORDINATE, Z 

Figure 21 .- Discretization of air filled wave guide 

WITH ABSORBERS MOUNTED ALONG BOTH THE UPPER AND LOWER 
WALLS. 
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Figure 22.- Effect of upper and lower wall absorbers 

ON THE MAGNITUDE OF THE AXIAL FLUX OF ENERGY (POYNTING 
VECTOR) FOR A FIVE-MODE MODAL EXPANSION IN THE EN- 
TRANCE AND EXIT DUCTS, 




Figure 23.- Contour plots of magnetic field amplitude. 










Figure 2‘j.- Zoom capabilities of finite element code. 















Figure 25.- Multiple thin film wall absorber configuration (not 

TO SCALE). 
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Figure 26.- Effect of multiple thin film wall absorber 

ON THE MAGNITUDE OF THE AXIAL FLUX OF ENERGY (POYNTING 
VECTOR). 
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